THE GALILEO SATELLITE CONSTELLATION AND MODIFICATIONS TO 
THE INVERSE-SQUARE LAW FOR GRAVITATION 



J. Paramos^ and O. Bertolami^* 



o 
o 

(N 

o 
O 

oo 



o 

I 



> 

m 



o 



X 



Instituto de Plasmas e Fusdo Nuclear, 
Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 
^ e-mail address: paramos@ist.edu 
^ e-mail address: orfeu@cosmos.ist.utl.pt 



ABSTRACT 



We consider the impact of a power-law correction to 
the Newtonian potential, inspired by ungravity or 
f{R) extensions of the Standard Model, and draw 
conclusions on the possibility of measuring the rele- 
vant parameters through observables made available 
by the Galileo satellite positioning system. 

Key words: Power-law additions, Ungravity, f{R) 
theories. 



INTRODUCTION 



The Gahleo positioning system poses a great oppor- 
tunity, not only for the improvement and develop- 
ment of new applications in navigation monitoring 
and related topics, but also possibly for fundamen- 
tal research in physics. Indeed, together with the 
already deployed Global Positioning System (GPS) 
and Glonass, satellite navigation may be considered 
the first practical application where relativistic ef- 
fects are taken into account, not from the usual ex- 
perimental scientific point of view, but as a regular 
engineering constraint on the overall design require- 
ments. Indeed, effects arising from special and gen- 
eral relativity (GR) - gravitational blueshift, time 
dilation and Sagnac effect - may account to as much 
as ~ 40 fis/day, which is many orders of magni- 
tude above the accuracy of the onboard clock de- 
ployed in these systems. Moreover, the gravitational 
Doppler effect, of the order of Vn/c'^ ~ 10"^*^ (where 
Vat = GMe/Re is the Newtonian potential, G is 
Newton's constant. Me — 6.0 x 10^'* kg is the Earth's 
mass. Re — 6.4 x 10^ m is its radius and c is the speed 
of light) falls within the 1G~^^ frequency accuracy of 
current space-certified clocks, and must therefore be 
taken into account. In the GPS, this is done by im- 
posing an offset in the onboard clock frequency, while 



in the Galileo Navigation Satellite System (GNSS) 
this correction should be accounted by the receiver. 
For further details, the reader is directed to Refs. 
[1, 2, 3, 4] and references therein. 

This said, it is not clear as to what extent the ac- 
curacy of the Galileo positioning system may be im- 
proved - which is designed to offer pinpoint localiza- 
tion within an error margin of 1 to, against the 10 m 
margin of previous the GPS - so to provide clues 
to the nature of models beyond the current gravita- 
tional physics knowledge (see Refs. [5, 6] for updated 
surveys). In this study, we aim to establish bounds 
on the detectability of extensions to the recently dis- 
cussed unparticle extension of the Standard Model 
of particle interactions [7] that manifest themselves 
through power-law corrections to the Newtonian po- 
tential [8]. These corrections are associated to the 
exchange of modes that do not correspond to par- 
ticle states, usually referred to as unparticles. The 
exchange of spin-2 states lead to a correction to New- 
ton's potential usually dubbed as ungravity. 

This paper is organized as follows: firstly, we assess 
the main relativistic effects that are present in the 
GNSS. We proceed and consider the possibihty of 
measuring the discussed corrections using the GNSS. 



MAIN RELATIVISTIC EFFECTS 
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2.1. Frame of reference 



Assuming that all time-dependent effects are of cos- 
mological origin, and hence of order Hq^, where 
Hq is Hubble's constant, one may discard these as 
too small within the timeframe of interest; hence, 
one assumes a static, spherically symmetric scenario, 
posited by the standard Schwarzschild metric. In 
isotropic form, this is given by the line element 
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where dV = dr"^ + dJl^ is the volume element, and 
V is the gravitational potential. In the standard GR 
scenario, the latter coincides with the Newtonian po- 
tential V = Vn = -{l-\- T.-2=iJn)GME/r, where the 
Jn multipolcs account for the effect of geodesic per- 
turbations and density profiles. 

To this picture, one must introduce the rotation of 
the Earth with respect to this fixed-axis reference 
frame, with angular velocity w = 7.29 x 10~^ rad/s; 
by doing a coordinate shift t' — t, r' = r, 9' = and 
(j)' = (j) — ujf' , one gets the Langevin metric, given by 
the line element 



Finally, if one gets back to a non-rotating frame, the 
metric is given by the line element 
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2.2. Constant and periodic clock deviation 

One may now consider the difference between the 
time elapsed on the ground and the satellite clock; 
keeping only terms of order c~^, one finds that the 
proper time increment on the moving clock is given 
approximately by 
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where, for simplicity, primes were dropped. Asides 
from a non-diagonal element, one obtains an addi- 
tion to the gravitational potential, which could be 
viewed as a centrifugal contribution due to the rota- 
tion of the reference frame. One can then define an 
effective potential $ = 2V — (ujrsinO)'^; the parame- 
terization of the Earth's geoid is obtained by taking 
the multipole expansion of V up to the desired order 
and finding the equipotential lines $ = $0 (the lat- 
ter being the value of $ at the Equator), and solving 
for r{6, (j)). 

In the above line elements, the coordinate time coin- 
cides with the proper time of an observer at infinity. 
However, since one wishes to evaluate the ground 

to orbit clock synchronization, it is advantageous to 
rewrite the metric in terms of a rescaled time coordi- 
nate, which coincides with the proper time of clocks 
at rest on the surface of the Earth; this is best imple- 
mented by resorting to the above-mentioned geoid, 
since its definition as an equipotential surface <I> = $0 
indicates that all clocks at rest with respect to it will 
beat at the same rate; hence, rescaling the time co- 
ordinate according to f — > (1 -|- $o/c^)i) one gets the 
metric given by the line element 
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Considering an elliptic orbit with semi-major axis a, 
and taking V = Vn — GMe/t, this may be recast 
into the form [1] 



dr = ds/c = 
3GMe 
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The constant correction terms in the above amount 

for the GNSS, and -4.4647 x 10-^°, for the GPS; this 
indicates that the orbiting clock is beating faster, by 
about 41 fis/day, for the GNSS, and 39 jis/day, for 
the GPS. For this reason, the GPS has a built in fre- 
quency offset of this magnitude, while the increased 
computational capabilities made available to current 
and future receivers of the GNSS leave this correc- 
tion to the user. The residual periodic corrections, 
proportional to 1/r — 1/a, have an amplitude of or- 
der 49 ns/day, for the GNSS, and 46 ns/day, for the 
GPS. 



2.3. Shapiro time delay and the Sagnac ef- 
fect 



The so-called Shapiro time delay, a second order rel- 
ativistic effect due to the signal propagation is given 

by[i] 
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where we have integrated over a straight line path 
of (proper) length I. Evaluating this delay, one con- 
cludes that this effect amounts to 6.67 x 10~^^ s. 

Also, one must consider the so-called Sagnac effect, 
which arises from the difference between the grav- 
itational potential V and the effective potential 
when proceeding from a non-rotational to a rota- 
tional frame. Hence, one gets the additional time 
delay 



^tsagnac = " / r"^ d(j) = — dA^ , (9) 

J path ^ J path 

where dA^ is the ortho-equatorial projection of the 
area element swept by a vector from the rotation 
axis to the satellite. For the GNSS, this yields a 
maximum value of 153 ns, while for the GPS, one 
gets 133 ns. 

One conchides this section by repealling the main ef- 
fects affecting the considercid global positioning sys- 
tems: a frequency shift of order 10^^'' and a prop- 
agation time delay (Shapiro plus Sagnac effect) of 
the order 10~^ s. In what follows, one will com- 
pute the additional frequency shift and propagation 
time delay induced by ungravity inspired corrections 
to the Newtonian potential, and compare the results 
with the above quantities, plus the frequency accu- 
racy of 10"^^ and the time accuracy of Galileo, of 
order 10~^ s, which corresponds to a optimistic spa- 
tial accuracy of 30 cm. 



more usual Yukawa parameterization of modifica- 
tions to gravity, thus enabling one to probe poten- 
tial candidates for extensions and modifications of 
GR. Furthermore, some current proposals give rise 
to the existence of induced power-law effects at as- 
trophysical scales. As mentioned, corrections of this 
type arise from the exchange of spin-2 unparticles 
from which one can set bounds, for instance, from 
stellar stability considerations [13] and cosmological 
nucleosynthesis [14]; in this context, the exponent is 
related to the scaling dimension of the unparticle op- 
erators du through n = 2du — 2, and the lengthscale 
R reflects the energy scale of the unparticle interac- 
tions, the mass of exchange particles and the type of 
propagator involved. 

Another relevant candidate for detection through 

a power-law parameterization is the so-called f{R) 
class of theories [15], which relies upon a generaliza- 
tion of the Einstein-Hilbert action (including a non- 
minimal coupling of geometry with matter [16]); in 
an astrophysical setting, these give rise to power-law 
additions to the Newtonian potential, which may be 
relevant in the context of the galaxy rotation puzzle 
[17, 18]. 

In the above Eq. (10), the Newtonian potential $jv 

is recovered by setting i? = (for positive n) or 
i? — > (X) (for negative n). The limit n — > is ill- 
defined, since then the additional term Vp does not 
vanish, but is identified with the Newtonian poten- 
tial, Vp = ^n: this indicates that one should rewrite 
the gravitational constant in terms of an effective 
coupling, leading to the full potential 



3. DETECTION OF A POWER-LAW AD- 
DITION TO THE NEWTONIAN PO- 
TENTIAL 



Having studied the measurability of an array of hy- 
pothetical perturbations to the Newtonian potential 
in a previous work [9] (including the effect of the 
cosmological constant, of a Yukawa term [10] and a 
constant acceleration such as the one reported in the 
so-called Pioneer anomaly — see Ref. [11] and refer- 
ences therein for a review, and Ref. [12] for a conven- 
tional explanation), one now considers a power-law 
perturbation of the form 



Vp = 



GMe 



(10) 



where n is a (possibly non-integer) exponent and R 
is a characteristic length scale arising from the un- 
derlying physical theory. 

A power-law perturbation to the Newtonian poses an 

interesting scenario for several reasons: phenomeno- 
logically, it provides a convenient alternative to the 
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with 



Gp = 



a 



(12) 



where i?o signals the distance at which the full 
gravitational potential matches the Newtonian one, 
$(i?o) = ^w(i?o)- 

This additional length scale Rq should arise as an 
integration constant when solving the full field equa- 
tions derived from the fundamental theory of grav- 
ity beyond the considered power-law correction. For 
simplicity, this study assumes that (R/Rq)" <§; 1, so 
that this may be safely discarded — at the cost of 
neglecting the regime n ^ 0, for which this condi- 
tion fails (notice that this is complementary to the 
approach followed in a previous study [13]). Hence, 
the siibsequent section uses Gp = G freely, albeit 
with the above caveat in mind. 
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3.1. Relative frequency shift 

The relative frequency shift of a signal emitted at a 
distance from the origin r = Re + h (for the GNSS, 
h = 23.222 X 10^ km) and received at a distance 
r = Re is given by 
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Hence, one may compute the additional frequency 
shift induced by the extra contribution to the poten- 
tial, through 
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where a dimensionless ratio ^ = R/ Re has been de- 
fined. One clearly has two asymptotic regimes: for 
n » — 1, one may approximate the above by 
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Figure 1. Contour plot for the relative frequency de- 
viation e as a function of^ = R/Re and n, with con- 
tours for e = 10~^^ (solid line), 10~^^ (long dash) 
and 10~^^ (short dash), and allowed region grayed 
out. 



3.2. Propagation time delay 

Besides the relative frequency shift discussed above, 
a power-law correction to the gravitational Newto- 
nian potential also induces a modification to the 
propagational time delay, given by 



: 6.96 X 10"^^" 



(15) 



so that, considering the accuracy e/^ = 10^^^ of the 
Galileo constellation, one has ^" < 1.44 x 10~^. Con- 
versely, if n -C — 1, one obtains 
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^ ) = -1.50 X 10-^(0.220" , (16) 
J J p 



Equating this to the 10 ^ s time resolution of the 
GNSS, yields the constraint 



which, upon comparison with the accuracy e = 
10-12, yields (0.22^-" > 1.50 x 10^ Since the r.h.s. 
is larger than unity, this imposes the rather strong 
bound ^ > (0.22)-i ~ 4.64 when n -co. 

In the vicinity of n = —1, one expands Eq. (14) as 



=1.07x10-^^" + ' 



/ 



(17) 



and the considered accuracy yields ^ > 1.07x 10^|n-|- 
1|. The different regimes are depicted in Fig. 1. 



1 - 0.22" 



n 



< 7.52 X 10 



-16 



(19) 



Analogously to the previous discussion, for the ex- 
treme negative power-law regime n <C this expres- 
sion reads |0.22C|" < 7.52 x 10"i*'|n|. Since, for 
sufficiently large (negative) n, the r.h.s. is larger 
than unity, one recovers the stronger bound, ^ > 
(0.22)-i ~ 4.64, obtained above. 

The allowed region for the ^, n parameters is de- 
picted in Fig. 2. 
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Figure 2. Contour plot for the propagation time delay 
Afp as a function ofS, = R/Re andn, with contours 
for Atp = 10"^ s (solid line), 10~^^ s (long dash) 
and 10~^^ s (short dash), and allowed region grayed 
out. 

4. CONCLUSIONS 



In this contribution, one has assessed the possibil- 
ity of detecting signals of new physics through the 
use of the GNSS. This application could be valuable, 
as any unexpected new phenomenology could pro- 
vide further insight into what lies beyond the Stan- 
dard Model of particle interactions and GR. One has 
specifically looked at the propagation time delay and 
frequency shift induced by a putative correction to 
the Newtonian potential, namely one of the a power- 
law type. This study complements the one where the 
GNSS was used to examine a constant cioiitribution, 
cosmological constant induced and Yukawa-type cor- 
rections to the Newtonian potential [9]. 

As can be seen from Figs. 1 and 2, one cannot lift the 

degeneracy in the {RE,n) parameter space through 
an intersection of the relative frequency shift and 
propagation time delay observables, since these pro- 
duce rather similar allowed regions for the mentioned 
parameters. 

Asides from the obtained range for these quantities, 
the main result of this study lies in the exclusion of 
a definite range of values for R, namely those lying 
between < i? < 4.64i?£:. Recall, however, that 
condition {R/RqY < 1 kads to the limit Rq < R 
when n — > — (X), yielding an additional requirement 
for the normalization constant Rq in the case of a 
large negative exponent n. Also, the above results 
are not valid in the vicinity of n = 0, due to the lack 
of knowledge concerning the effective gravitational 
coupling Gp. 
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